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I. 



INTRODUCTION 



Every product eventually fails and, sooner or later, is 
no longer usable for its original purpose. The rate at which 
this deterioration process proceeds is called the failure rate. 
One is interested in predicting the time to failure for a 
component or an entire system. It is impossible to predict 
the time to failure precisely, but one can model this process, 
viewing the time to failure as a random variable which has 
a probability distribution that can be specified. There 
exist numerous life distributions which provide realistic 
probabilistic descriptions of the behavior of components or 
whole systems. 

This paper is concerned only with the exponential distri- 
bution which plays a fundamental role in reliability. This 
distribution has many properties which facilitate the modelling 
process, and so encourage its application in many fields. 

Where not directly applicable, it is often used as a bound. 

According to the exponential distribution, a component 
has a constant failure rate. An equivalent property is that 
the remaining life of the component is independent of its age: 
the "memoryless" or "no aging" or "no wear" property. This 
means that a used component is always as good as a new one. 

Convolving life distributions for components to get the 
final life distribution for the system is rather difficult. 
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even for the exponential distribution. For an increasing 
number of components it's no longer possible to compute a 
closed form life distribution. This paper will show how to 
use a reliability shorthand notation to get analytical results 
for small basic systems, and later use an algorithm to com- 
pute numerical results with an Apple 11+ microcomputer. The 
Pascal programs are provided in the appendices. 
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II. RELIABILITY SHORTHAND 



A. NOTATION AND DEFINITIONS 

The shorthand notation "EXPtAjJ + EXP(A 2 ]" denotes the 
distribution for a random variable T^+T 2 . T^ has an exponen- 
tial distribution with failure rate A^ and T 2 has an exponen- 
tial distribution with failure rate A 2 * T-^ and T 2 are assumed 

to be independent . 

In the following T always denotes a random life length 

and 



F ( t ) = Probability (T > t) = 1 - F(t) 

is the survival function as a function of time. 

The convolution of life distributions corresponds to 
the summation of components with independent random lives 
to get the life distribution of a system. 
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III. RELIABILITY FOR A ONE COMPONENT SYSTEM WITH SPARES 



A. GENERAL DESCRIPTION 

There are two ways for a simple system, consisting of 
one active component and one spare, to survive a specified 
time t. Either the active component survives time t, in 
which case the spare need not be used, or the active com- 
ponent fails at some time s (0 <_ s <_ t) and the spare takes 
over and doesn't fail in the interval from s to t and thus 
completes the mission. 

Let F^(t) and F2(t) be the exponential survival functions 
for component one and two respectively, and f^(t) and f2(t) 
their respective densities. Then the survival function F(t) 
for the system is denoted by: 

t 

F (t ) = F (t) + / F- (t-s) f. (s)ds 

1 0 

where the integrand shows the probability that component one 
fails at time s (f^(s)) and component two takes over for 
at least the length of time t-s. The following shows some 
examples how to get the survival function in closed form for 
different basic convolutions. 
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B. BASIC SURVIVAL FUNCTIONS 



1 . System with N Identical Components 

The shorthand notation for this system is: 

EXP [A] + EXP [A] + ... + EXP [AJ 

Whenever an active component fails a spare takes over, as 
long as there is a spare. When the n-th component fails, the 
whole system stops functioning. 



For n = 2 : 



t 

F(t) = F x (t) + / F 2 (t-s)f 1 (s)ds 




0 



( 1 + A t) e 



-At 



t > 0. 



The density function for this survival function is: 



f (t) 



dF x (t) 



dt 



d [ ( 1+At) e At ] 
dt 



-[ (l+At) (-Ae At ) + e Xt A] 
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For n = 3: 



The survival function for the two component system 
will be convolved with the survival function of the next 
component to obtain the survival function for the three 
component system. 

F(t) = (l+At)e' Xt + / e" X (t " s) A 2 se" As ds 

0 

which leads to: 



F ( t ) 



( 1 + A t + 



(At) 2 , .-At 
2~l e ' 



t > 0 . 



For all n: 

Continuing in this manner one gets the general solution 
for the survival function of this kind of system: 



F ( t ) 



n . 

I (At) 1 " 1 
i=l 



( i-1) 




t > 0 . 



This class of convolutions always gives the well-known 
Erlang [n, A] survival function. 

2 . System with N Dissimilar Components 

The shorthand notation for this system is denoted 

by: 



EXPUJ + EXP[A 2 ] + ... + EXP[A r ] 
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For n 



2 : 



The survival function for the active component is 



-A, t 

F x (t) = e 1 , 



t > 0. 



The survival function for the spare is 



-A-t 

F 2 (t) = e z , 



t > 0 



The survival function for this system is 



F(t) = F (t) + / F- ( t-s) f . (s) ds 

i o 1 



-A..t t-A„(t-s) -A-s 

= e +/e A,e J ‘ds 



-A,t A. -A t t -(A -A_)s 

= e 1 + t — =-? — e z f ( A . - A _ ) e 1 z ds 
A 1 “ A 2 0 



-A t " A 2 t /1 (X 1 X 2 } 

= e + ? r — e (1 - e ) 

A 1 A 2 

-A . t -At 

(A 1 -A 2 )e + A x e 2 - A.^ x ) 

( Af “ X 2 > 

-A-t - A 1 1 

(A x e - A 2 e ) 

( A]_ “ A 2 ^ 



t > 0 . 
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For n 



= 3 : 

The resulting survival function for the two component 
system will be convolved with the survival function of the 
next component to obtain the survival function for the three 
component system. 

The result is 



F(t) = 



A 2 A 3 



( A 3 - A j_) ( A 2 — A 1 ^ 



-A x t 
e + 



A 1 A 3 



A 2 t 



( A 3 -A 2 )(A r A 2 ) 



A 1 A 2 



- A 3 t 



(A 3" A 2 } (A 3 _A 1 ) 



For all n: 

Continuing in this manner one gets the general solution 
for the survival function of this kind of system: 



F (t ) 



n 



1 

i=l 



3 

n (A. -A.) 
jri 3 1 



e 



-A . t 
1 



The two general survival functions for systems built from 
components which all have dissimilar failure rates or all 
have the same failure rate are very useful for computations 
of survival functions of systems which are formed from sub- 
systems of the above two kinds. 

It will be shown in the following section how to use 
the above results to get a closed form for the survival 
function of a system which can not be computed with only one 
of those general formulas . 
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3 . Systems With Identical and Dissimilar Components 



Shorthand notation: 

EXP[A 1 ] + EXP[A 1 ] + EXP[A 2 ] + EXP[A 2 ] 

Let the random variables T be distributed in the following 
way : 



T x ~ EXPUJ, T 2 ~ EXP[X 1 ], 

T 3 - EXP[A 2 ], T 4 - EXP[A 2 ]. 

There are different approaches available for computing the 
survival function for this system. 

One can compute the survival function for the first 
two and then the second two components, using the general 
formula of Section III.B.l and then convolve those two sur- 
vival functions; or one can compute the survival function of 
the first two components and then add the other components 
one at a time. The latter approach will be used here to 
give a general idea of the procedure, and to lead to an 
algorithm for the numerical evaluation of the reliability of 
systems consisting of so many components that a closed form 
is no longer attainable or not attainable in an efficient 
manner. 

The first two components (T^ = T^ , T 2 = T^ + T 2 ) 
have as shown in III.B.l the survival function 
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t 



F , (t) = F T ,(t) 


L- 

+ / 


F (t-s) f T , ( s) ds 


2 L 1 


0 


i 2 T 1 






•A, t 


= (a 11 + 


a 12 t)e 


h 

ft 

1 v 

O 


with : 






a n - 1 


and 


a 12 “ X V 



Adding the next component leads to T' = T. + T_ + T- 

J X Z j 

and the survival function 



F T'(t) 



Fm.(t) + / F , (t-s) f , (s) ds 

2 0 1 3 l 2 



-A,t t -A^(t-s) m — A s 

= (1 + A 1 t)e + / e 

0 



2 A?se ds 



which leads after the integration to 



F T , (t) 
3 



( A 2 ~ 2A 1 A 2 ) + (A A 2 - A^A 0 )t -A^t 

- e + 

( a 1 ~a 2 ) 



-A 2 t 



(A i-a 2 ) 



or in the form of equation (1) to 



with 



-A,t ~ A - 1 

F T ,(t) = (a i;L +a 12 t)e + a 21 e 



l ll 



A 2 ~ 2 A]_A 2 



(Ai~A 2 ) 



2 ' 12 



(X 1 ~X 2 ) 



21 
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+ T + and 



Adding the last component leads to = 
the survival function: 



F t , (t) 

L 4 



F T .(t) + / F , (t-s)f ,(s)ds 

3 0 l 4 l 3 



-A.t -At 

= (a u+ a 12 t)e + a 2] _e 



t -A 2 (t-s) -As 

+ J e [(a 12 -a n X 1 -a 12 X l s)e 



-A 2 s 

- ^2 a 21 e ] ds 



Integrating and collecting terms gives the final survival 
function for the whole system: 



-A, t -At 

F T ^(t) = (a i;L +a 12 t)e + ( a 21 +a 22 t) e 



t > 0 



with 



l ll 



X 2- 3X 1 X 2 
( a 2 ~ ^i') 



A A 2 

A 1 A 2 



12 



(a 2 -a 1 ) 



2 ' 



21 



3 2 

X 1 " 3X 1 X 2 

( A 2 — A x ) 3 



X 1 X 2 
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(a 1 -a 2 ) 



This section demonstrates that one can get the survival 
function by adding one random variable at a time and thus 
computing the new survival function. Each summation then 
gives a new or updated a. .. Even though the integration is 
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not shown explicitly, one can imagine that this method of 
convolving survival functions is very tedious and it is very 
easy to make mistakes. So the next section introduces an 
algorithm to compute all the a^j , which then makes it easy 
to evaluate the survival function at any time. 



C. ALGORITHM TO COMPUTE THE CONVOLUTION OF EXPONENTIAL LIVES 
The shorthand notation to be solved is 



EXP [ A x ] + EXPUJ + . 

+ exp[a 2 ] + expu 2 ] + . 



. + EXPU-J 

. + exp[a 2 ] 



(n.^ terms) 
(n 2 terms) 



+ EXP t A k ] + EXP[A k ] + 



. + EXP[A k ] 



(n k terms) 



The general solution for the survival function is of the form 



— A. -| t ~ A. t 

F(t) = A 1 (t)e 1 + A 2 (t)e + ... + A k (t)e , t > 0 



with 



2 


^3 


(ni-1) 


11 + a 12 t + a 1 1 6 


+ a i4 t + 


••• + a ln 1 t 


2 


^3 


( n 2 -1) 


21 + a 22 t + a 23 t 


+ a 24 t + 


••• + a 2n 2 t 



A k (t) ■ a kl + a k2 t + a k3 t2 + a k4 t3 + ••• + a kn k t 



(n k -l) 



An example shows how to get the survival function for the 
following shorthand notation: 
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EXPUJ + EXP[X 1 ] + EXPU^ 
+ EXP[X ] + EXP[X 2 ] 

+ EXPU 3 ] 

Survival Function; 



F ( t) = 



2 1 ^" 

(a i:L +a 12 t+a 13 t )e + 



-X-t -X-t 

^ a 21 +a 22 t ^ e + a 31 e 



t > 0 . 



An algorithm to compute all the factors a^ , a^ 2 , . .., a^ 
and thus make it possible to convolve all the different 
exponential functions was earlier introduced by Gursel [Ref. 
1] , and is given here in a slightly modified form. This 
algorithm adds one exponential random variable in each run. 
The above example can be solved in this way by looping six 
times through this algorithm. 

Notation for the algorithm; 



K = number of dissimilar failure rates (3 in above 
example) ; 

1^ = i-th type failure rate, i = l,...,k; 

nn . = number of random variables having the 9-th 

type failure rate; 



1. = failure rate of the ie-th lambda, ie = 1,...,K 

re 

n. = current number of random variables having the 
1 i-th failure rate, n^ = l,...,nn^; 



a^j = j-th coefficient of the i-th polynomial 

ie = index for current polynomial and currently 
entering life. 
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Algorithm to compute the coefficients 



Initialization : 


a .. 
Dk 


= 0, for all j,k where j 
k = 1 , . . . , nn . 

i 


1 , . . . , K , 




n. 

i 


= 0, for all i where i = 


• 

• 

• 

l — 1 


Input : 


X . 
i 


for all i where i = 1,... 


,K 




nn . 
l 


for all i where i = 1,... 


,K. 


After the first 


run : 


a-^^ — 1, n^ — 1 





n 



ie 



n. +1, until n. 
le le 



nn 



ie 



1. Update the coefficients a. , for k = 2,...,n. 

1 e r k ig 



a. . = A. *a. -./(in-l) , where m = n. -j for 

ie,k le ie,ra-l / le J 

j = 0 , . . . ,n. -1 
J ' le 



2. Update the coefficients a. , 

16 / i 



a . -« = a . 

ie,l le 



> ( j -11 iA ie a ii 



n . , n . 

l a . n A . l 

1 + I -A” T + l 3 

' i^ie U i v ie } j=2 (A -A. ) 3 

i=l i^ie 1 ie 



3. Update the other coefficients a. . for all i,k where 

i <> ie , n^ > 0 

for i = 1 , . . . ,K, 
k = 1 , . . . ,n^ 



a . = a . A. /(A. -A.) for all i ^ ie, 

i,n. i,n. ie le i 

i i 



i,n. 



(a . A . -ma . ) / ( A . -A . ) 

i,n. ie l ,m ' ' ie i 

' l 



for all i where 
i ^ ie and n. > 1 
for m = n . - j , 
j = 1, . . . , (n^-1) 
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D. PROGRAM FOR NUMERICAL EVALUATION OF THE RELIABILITY 
FOR A ONE COMPONENT SYSTEM WITH SPARES 

The above algorithm is used to compute the reliability of 
a one component system with similar and/or dissimilar spares 
for any given time. The input to the program, which is 
written in PASCAL and listed in Appendix A, will be explained 
here with the same example as above: 

EXP[X X ] + EXP [ A x ] + EXPU^ 

+ expu 2 ] + exp[a 2 ] 

+ exp[a 3 ] 

Let's assume that A^ = 0.1, A 2 = 0.2, A 3 = 0.3 and we want 
to know the reliability for t = 0, 10, 20, ..., 90 hours. The 

prompts of the program are typed in large letters . 

"HOW MANY DISSIMILAR LAMBDAS DO YOU WANT TO USE?" 

3 

"WRITE DOWN ALL THE DISSIMILAR LAMBDAS WITH AT LEAST ONE 
BLANK IN BETWEEN (INPUT NUMBERS < 1 WITH A LEADING ZERO) : " 

0.1 0.2 0.3 

"WRITE DOWN THE NUMBER OF EACH OF THE LAMBDAS WITH AT LEAST 
ONE BLANK IN BETWEEN:" 

3 2 1 

"HOW MANY TIME POINTS DO YOU WANT TO BE COMPUTED? (<=20) :" 

10 

"AT WHAT TIMES DO YOU WANT THE SURVIVAL FUNCTION TO BE 
EVALUATED? (TIME SHOULD HAVE THE SAME DIMENSION AS THE 
LAMBDAS AND SHOULD BE AT MOST 20 INTEGERS BETWEEN 0 AND 70. 
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THE UPPER LIMIT 70 IS REQUIRED ONLY, IF A GRAPHICAL 
PRESENTATION IS DESIRED LATER) 

0 10 20 30 40 50 60 70 80 90 

"DO YOU WANT A SIMULATION TO BE DONE PARALLEL TO THE 
COMPUTATIONS? (SERIOUS SIMULATION NEEDS A LOT OF TIME!) HIT 

1 FOR "YES" AND "0" FOR NO:" 

1 

"HOW MANY NUMBERS DO YOU WANT TO USE AS BASIS FOR THE 
SIMULATION? :" 

1000 



After the computations the program will produce the following 
output : 

FOLLOWING LAMBDAS WERE ENTERED: 0.1 0.2 0.3 

# OF EACH LAMBDA: 321 



Results : 



RELIABILITY 


( 


T 


> 


0 


) 


= 


1.000000 


(SIM 


= 


1.000) 


RELIABILITY 


( 


T 


> 


10 


) 


=: 


0.995799 


(SIM 


= 


0.996) 


RELIABILITY 


( 


T 


> 


20 


) 


= 


0.925259 


(SIM 


= 


0.925) 


RELIABILITY 


( 


T 


> 


30 


) 


= 


0 . 739626 


(SIM 


= 


0.728) 


RELIABILITY 


( 


T 


> 


40 


) 


=: 


0.508910 


(SIM 


= 


0.489) 


RELIABILITY 


( 


T 


> 


50 


) 




0.311271 


(SIM 


= 


0.334) 


RELIABILITY 


( 


T 


> 


60 


) 


= 


0.174439 


(SIM 


= 


0.173) 


RELIABILITY 


( 


T 


> 


70 


) 


= 


0.091597 


(SIM 


= 


0.068) 


RELIABILITY 


( 


T 


> 


80 


) 




0.045784 


(SIM 


= 


0.036) 


RELIABILITY 


( 


T 


> 


90 


) 


= 


0.022028 


(SIM 


= 


0.020) 



"DO YOU ALSO WANT A GRPAHICAL PRESENTATION OF THE DATA? (HIT 
"1" FOR "YES" AND "0" FOR "NO"):" 
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1 



The graphical output will not be shown here, but will be 
generated from the program if desired. 

Some notes on the program: 

1. The program is set up to allow 20 dissimilar lambdas. 

2. It is necessary to type leading zeros in PASCAL if one 
inputs decimal numbers. 

3. The simulation generates n numbers (in this example 
n = 1000) for every random variable, to produce n 
observations on the whole system. 

4. Because of 

the behavior of the exponential function, 
the sequence of the lambdas in the input, 
the method of adding always only one random variable 
at a time, 

the possible closeness of the values of the lambdas 
and 

the possibly large number of the dissimilar lambdas, 
the program might have to compute with very large and 
small numbers and could produce erroneous results. 
Therefore a simulation is provided within the program 
to see if the computations are close to the simulation. 
Fixed limits which avoid errors can't be given, because 
all the possible errors are closely related. Fewer 
dissimilar lambdas for example would allow that their 
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values could be closer together, or that less care is 
needed with the sequence in the input, etc. 

But after a long test phase and a comparison with results, 
produced with the IBM 3033, some hints can be given: 

1. The input sequence of the lambdas should go from the 
smallest to the largest, 

2. There should be no more than 20 lambdas altogether, 

3. The difference in value of the dissimilar lambdas should 
be larger than 0.05. 
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IV. RELIABILITY FOR A SYSTEM OF ACTIVE COMPONENTS AND SPARES 



A. MIXTURE OF CONVOLUTIONS 

Before solving systems with active components and spares, 
the mixture of convolutions will be explained briefly. Let's 
assume that we want to compute the reliability of a parallel 
system with two different components, like that in Figure 1. 




Figure 1. Simple Parallel System 

The system starts with a life distributed as EXPIX^^]* Then 
there are two different possibilities: either component #1 

fails with probability = X- l /(A ;l +X 2 ) and the residual life 
is distributed as EXP[A 2 ]. This is the upper branch of 
Figure 2. Or component #2 fails with probability 
P 2 = X 2 // ^l + ^2^ and the residual life is distributed as 
EXP [X ± ] . 
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Figure 2. Shorthand Notation for Parallel System 



Using the "MIX"-notation, the life can be written in the 
shorthand notation as: 



EXP[A 1 +A 2 ] + MIX{P 1 *EXP [A 2 ] ,P 2 *EXP [A-J } . 

Because of the distributive law, which is valid for these 
mixtures, one can transform this equation to 

MIX{P 1 *EXP [A 2 ] + EXP [ A jl+A 2 ] ) ,P 2 * (EXP [A,] + EXPlAj+Aj])} 

The first and the second expressions in parentheses multiplied 
by their respective probabilities, is a mixture of convolu- 
tions. The following section demonstrates how those mixtures 
can be solved numerically even if an analytical solution is 
not feasible. 
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B. DESCRIPTION OF THE SYSTEM 



This section is concerned with computing the reliability 
of a system of components. The different components might 
be supported by spares, which have the same or a different 
failure rate as the original component. The following example 
will be used to illustrate the idea of describing the system 
in the reliability shorthand notation. The listing of a 
PASCAL program for these types of systems is provided in 
Appendix B. The program incorporates the method of program 
#1 to add the components, one at a time, to get first the 
convolution of the exponential distributions and then via 
the mixture the reliability of the whole system. 

Example : 

The system consists of two parallel components (#1 and 
#2) in series with a third component (#3) . A spare is pro- 
vided for component #1. The components have the failure rate 
A^, A^ and A^ respectively. The spare has the same failure 
rate as the component it can replace. 




Figure 3. Combined Serial/Parallel System with Spare 
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C. SHORTHAND BRANCHING DIAGRAM 



A l +A 2 +A 3 / 

X 

A 2 




EXP[A 2 +A 3 J (1) 

EXP[A 2 +A 3 J (2) 

ZERO (3) 

EXPIAj+A^ (4) 

ZERO (5) 



^ZERO (6) 



Figure 4. Shorthand Branching Diagram for Serial/ 
Parallel System with Spare 



The numbers in the numbered parentheses in the figure 
above describe the different endpoints (leaves) in this nota- 
tion. The connections between the original system (root) 
and the leaves will be called branches. In the following: 
the connection from the root to leaf #1 is called branch #1, 
from the root to leaf #2 branch #2, etc. 

The original system has a life distributed as EXP lA^+A^A^J . 
Three different things could happen now. Either component 
#1 fails with probability A^/ (A-^+A^+A^ or component #2 fails 
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with probability A 2 / (A-^+A 2 +A 3 ) or component #3 fails with 
probability A 2 / ( A ^+A 2 +A 3 ) . Let's continue in branch #1 and 
assume component #1 failed with the above probability. Now 
the spare can replace the failed component and the system's 
life is again distributed as EXP [A -^+A 2 +A 3 ] . The system is now 
in the same state as at the beginning and the same arguments 
hold. Let's assume new component #1 fails with probability 
A^/ (A 1 +A 2 +A 3 ) , then the remaining life is distributed as 
EXP[A 2 +A 3 ]. If component #2 fails (the probability of this 
failure is A 2 (A-^+A 2 +A 3 ) ) , the remaining life is distributed 
as EXPtA-^+A^j. The third possibility, component #3 fails with 
probability A^/U-^+A^A^) , brings the failure of the whole 
system, denoted by the "ZERO "-notation. This says that there 
is a remaining life of zero. The reason why this third branch 
isn't simply omitted is that even if there is no remaining 
life, there is a specified probability that this event occurs. 
And this has to be taken into account to compute the total 
life of the system. 

The other three branches of the above figure are developed 
analogously . 



D. THE SURVIVAL FUNCTION 

Let's denote the survival function of branch #1 by F-^(t), 
and the probability that the events in branch #1 occur in 
this special sequence as P^. Then the life of branch #1 is 
distributed as P^* (EXP [A 1 +A 2 +A 3 ] + EXP [Aj+A^A^ + EXPIA 2 +A 3 J ) 
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and the survival function of this system is the following 
mixture : 

F(t) = MIX{p 1 *F 1 (t) »P 2 *F (t) ,p 3 *F 3 (t) ,p 4 *F 4 (t) ,p 5 *F 5 (t) ,p 6 *F g (t) } 

For the computation of the coefficients for every single 
survival function the algorithm of Section III.C can be used 
again . 



E. PROGRAM FOR THE NUMERICAL EVALUATION OF THE RELIABILITY 
FOR SYSTEMS WITH SEVERAL ACTIVE COMPONENTS AND SPARES 

This section describes the program for computing the 

reliability for systems as described above. While reading 

this part, one should look at Figure 4. The program is 

interactive. Only the input of Figure 4 will be explained 

here . 

Prompt Input 

"How many branches?" 6 CR (Carriage return) 

1. Branch : 



1. 


Element : 


2 


CR 


3 


CR 


0 


CR 










2. 


Element : 


1 


CR 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


3. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 






4. 


Element : 


1 


CR 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


5. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 






6 . 


Element : 


0 


CR 


















2. 


Branch : 






















1. 


Element : 


1 


CR 


3 


CR 


0 


CR 










2. 


Element : 


2 


CR 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


3. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 
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) 



4. 


Element : 


1 


CR 


1 


CR 


2 


CR 


3 


CR 0 CR 


5. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


6. 


Element : 


0 


CR 















3 . Branch : 



1 . 


Element 


• 

• 


0 


CR 


















2. 


Element 


: 


3 


CR 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


3. 


Element 


: 


1 


CR 


2 


CR 


3 


CR 


0 


CR 






4 . 


Element 


: 


1 


CR 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


5. 


Element 


: 


1 


CR 


2 


CR 


3 


CR 


0 


CR 






6 . 


Element 


; 


0 


CR 



















4 . Branch : 



1. 


Element : 


1 


CR 


3 


CR 


0 


CR 






2. 


Element : 


1 


CR 


1 


CR 


3 


CR 


0 


CR 


3. 


Element : 


1 


CR 


3 


CR 


0 


CR 






4 . 


Element : 


2 


CR 


1 


CR 


2 


CR 


3 


CR 0 CR 


5. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


6 . 


Element : 


0 


CR 















5 . Branch : 



1. 


Element : 


0 


CR 














2. 


Element : 


3 


CR 


1 


CR 


3 


CR 


0 


CR 


3. 


Element : 


1 


CR 


3 


CR 


0 


CR 






4. 


Element : 


2 


CR 


1 


CR 


2 


CR 


3 


CR 0 CR 


5. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


6. 


Element : 


0 


CR 















6 . Branch : 



1. 


Element : 


0 


CR 














2. 


Element : 


3 


CR 


1 


CR 


2 


CR 


3 


CR 0 CR 


3. 


Element : 


1 


CR 


2 


CR 


3 


CR 


0 


CR 


4. 


Element : 


0 


CR 
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The branches in the figure count from top to bottom and the 
elements in each branch from the leaf to the root. The first 
element in the first branch is EXP[A 2 +A 2 ] and one inputs only 
the two indices of the lambdas. The second element in the 
same branch is A^/ ( A 1 +A 2 +A 3 ) and one inputs first the index 
of the lambda in the numerator and then the indices of the 
denominator. After each element input a "0". If there is no 
sixth element as in branch #1 to branch #5, no fourth element 
as in branch # 6 , or no element as the ZERO elements in 
branches #3, #5, and # 6 , simply input a "0". Because of the 
data structure used in the PASCAL program, the system always 
has to be split up into pieces, so that there is only a single 
failure in each step. Because of the algorithm used, there 
shouldn’t be more than 20 elements in each branch to avoid 
errors . 
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V . SUMMARY 



The reliability shorthand is a convenient way to describe 
the life of different systems. This shorthand then is the 
underlying data structure to compute numerically the relia- 
bility for systems for which it is difficult or impossible 
to get the survival function in closed form. 
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APPENDIX A 



PROGRAM CONVOLUTION; 

USES TRANSCEND , APPLESTUFF ; 

CONST 

ASTERISKS *’ ; 

BLANK = ’ ’ ; 

DASH = 

LINE = ’ ! ’ ; 

TYPE 

ARR = PACKED ARRAY Cl.. 20 3 OF REAL; 

ARI = PACKED ARRAYC 1 . . 203 OF INTEGER; 

VAR 

A: PACKED ARRAY C 1 .. 20, 1 .. 203 OF REAL; 

L, PRO, REL, RELSIM : ARR; 

NINIT,NI,T: ARI; 

KK, JJ, IE, k , TMAX , SI MMAX , Y , X . COUNT , VAL: INTEGER; 

I, 1 1 , J , F ACT , NN , NKK , NN I , CHGR , CHS I M : INTEGER; 

I T , D .-INTEGER; (* INDEX FOR T-ARRAY *> 

SUM. RELI: REAL; 

PRINTVAR: CHAR; 

(* L = ARRAY OF DISSIMILAR FAILURE RATES *> 

(* NICI3 = NUMBER OF LIVES HAVING THE ITH FAILURE RATE *> 
(* AEI.J3 = COEFFICIENT FOR THE POLYNOMIAL *> 

(* T = ARRAY OF TIMES AT WHICH THE RELIABILITY WILL BE 
COMPUTED *) 

PROCEDURE ATTENTION; 

VAR 

DURATION, I; INTEGER; 

PITCH: ARRAYC 1 .. 103 OF INTEGER; 

BEGIN 

P I TCHC 13:= 15; 

PITCHC23 : = 17; 

PITCHC33 : = 19; 

PITCH! 43 : = 20; 

P I TCHC 53 : = 22; 

DURATION: = 30; 

FOR I : =1 TO 5 DO 

NOTE (PITCHC I 3, DURATION) ; 

END; 

PROCEDURE DECIDE2; 

BEGIN 

WRITELN; WRITELN; 

WRITELN ( ' DO YOU WANT TO CONTINUE WITH THE SAME PROGRAM? 
TYPE "U" : ' ) ; 

WRITELN (’DO YOU WANT TO COMPUTE MIXED CONVOLUTIONS? 

TYPE " XP2" :’); 

EXIT (PI ) ; 
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END; 

PROCEDURE INPUT; 

BEGIN 

WRITELN; WRITELN; 

GOTOXY (O, 12) ; 

WRITELN ('HOW MANY DISSIMILAR LAMBDAS DO YOU WANT TO USE? ’) 
READLN(K) ; 

WRITELN; 

WRITELN ('WRITE DOWN ALL THE DISSIMILAR LAMBDAS WITH AT 
LEAST ONE BLANK IN BETWEEN:'); 

WRITELN (’ (NUMBERS < 1 NEED A LEADING ZERO!)'); 

FOR I : =1 TO K DO 
READ (L C I 3 ) ; 

WRITELN; 

WRITE (’WRITE DOWN THE NUMBER OF EACH OF THE LAMBDAS WITH AT 
LEAST ONE BLANK IN’ ) ; 

WRITELN ( 'BETWEEN: ’ ) ; 

FOR I : = 1 TO K DO 
READ (NINITCI3) ; 

WRITELN; 

WRITE (’HOW MANY TIME POINTS DO YOU WANT TO BE COMPUTED 0 
( O 20 ) : ’ ) ; 

READLN (TMAX ) ; 

WRITELN (’AT WHAT TIME DO YOU WANT THE SURVIVAL FUNCTION TO 
BE EVALUATED?'); 

WRITELN( '( INTEGER TIMES BETWEEN 0-60 IN THE SAME UNITS AS 
THE LAMBDAS) ’ ) ; 

FOR I : =1 TO TMAX DO 
READ (TCI 3) ; 

WRITELN; 

WRITELN (’DO YOU WANT A SIMULATION TO BE DONE PARALLEL TO 
THE COMPUTATIONS?’); 

WRITELN (’ (SERIOUS SIMULATION NEEDS A LOT OF TIME!)’); 

WRITE (’HIT "1" FOR YES OR "O" FOR NO : ’); 

READLN (CHSIM) ; 

WRITELN; 

IF CHSIM = 1 THEN 
BEGIN 

WRITELN (’HOW MANY NUMBERS DO YOU WANT TO USE AS BASIS 
FOR THE SIMULATION?’); 

READLN (SIMMAX) ; 

END; 

WRITELN; WRITELN; 

END; 

PROCEDURE OUTPUT; 

BEGIN 

WRITE (’FOLLOWING ’,K:3,’ LAMBDAS WERE ENTERED :’); 

FOR I : =1 TO K DO 
WRITE(LCI3:5:2) ; 

WRITELN; 
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# OF EACH LAMBDA 



WRITE ( ' 

FOR I : =1 TO K DO 
WRITE <NINITCII:5) ; 
WRITELN; WRITELN; WRITELN; 
end; 



>; 



PROCEDURE INITIALIZE; 

BEGIN 
KK i=l; 

WHILE KK<=K DO 
BEGIN 

A IKK, II: =1.0; 

NI CKKI : = 1 ; 

J J : =1 ; 

WHILE J J < =K DO 
BEGIN 

IF JJOKK THEN 
BEGIN 

ACKK, II: =ACKK, 1I*LCJJI/ (L C J J I-LCKKI ) ; 
END; 

J J : =J J+l ; 

END; 

KK:=KK+1 ; 

END; 

KK: =K; 

IE: =1 ; 

END; 



PROCEDURE SUMMATION; 

BEGIN 

SUM: =0. O; 

FOR I : = 1 TO K DO 
BEGIN 

IF IOIE THEN 
BEGIN 

SUM: =SUM+AC 1 , 1 I*LC 1 1/ (LC I I-LCIEI > ; 

IF N I C 1 1 >=2 THEN 
BEGIN 

FACT : = 1 ; 

NKK: =NI C I I ; 

FOR II: =2 TO NKK DO 
BEGIN 

FACT : =FACT * ( 1 1 - 1 ) ; 

IF LCII>LCIEI THEN 

SUM: =SUM+AC I , I 1 1 *FACT*LC I El /EXP ( I I *LN ( 
LCII-LCIEI) ) 

ELSE 

BEGIN 

IF (II MOD 2) =0 THEN 

SUM:=SUM+AtI , I 1 1 *FACT*L C IEI /EXP ( 

1 1 *LN(ABS (LCII-LCIEI) > ) 

ELSE 
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SUM: =SUM+AC 1 , 1 1 3 *FACT *L C IE 3 / (-EXP< 
II*LN<ABS(LCI3-LCIE3> ) ) ) 

END; 

END; 

END; 

END; 

END; 

AC IE, 1 3: =SUM+AC IE, 13; 

END; 

PROCEDURE SUB1 ; 

BEGIN 
I : = 1 ; 

WHILE I< = K DO 
BEGIN 

IF IOIE THEN 
BEGIN 

ACI,NICI33:=ACI,NICI33*LCIE3/ (LCIE3-LCI3) ; 

IF NICI301 THEN 
BEGIN 
J : = 1 ; 

NNI : =NI C I 3-J ; 

WHILE J<= (NI C I 3-1 ) DO 
BEGIN 

AC I , NNI 3:=(LCIE3iACI, NNI 3 -NNI * AC I , NNI + 13)/ 
(LCIE3-LCI3) ; 

J:=J+1; 

NNI : =NI C I 3-J ; 

END; 

END; 

END: 

I : =1+1 ; 

END; 

END; 

PROCEDURE ONEATATIME; 

BEGIN 

IF NINITC IE3=NI C IE3 THEN 
IE: =IE+1 
ELSE 
BEGIN 

NI C IE3 : =NI C IE3 + 1 ; 

J : =0; 

REPEAT 

NN: =NI C IE3-J ; 

AC IE, NN3 : =LC IE3*ACIE, NN-1 3/ (NN-1) ; 

J : =J+1 ; 

UNTIL NN=2; 

SUMMATION; 

SUB1 ; 

END; 

END; 
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PROCEDURE RAND; 

VAR 

1 , 1 1 , 1 1 1 , COUNT , J : I NTE6ER; 

R, TSIM: REAL; 

BEGIN 

Ills =0; 

COUNT: =0; 

REPEAT 

TSIM: =0. 0; 

1 1 : = 1 ; 

FOR J : =1 TO K DO 
BEGIN 

FOR I : =1 TO NINITCJ] DO 
BEGIN 

R : =RANDOM/MAX I NT ; 

TSIM: =TSIM+ < — (LN (R) /LCII3) ) ; 

END; 

1 1 : = 1 1 + 1 ; 

END; 

IF TSIM>Tt IT 3 THEN 
COUNT :=C0UNT+1; 

III:=III+1; 

UNTIL 1 1 I=SIMMAX ; 

RELSIMC I T 3 : =COUNT/SIMMAX ; 

END; 

PROCEDURE SUB2 ; 

BEGIN 

RELI : =0. 0; 

IF TC IT 3=0 THEN 
RELI: =1.0 
ELSE 
BEGIN 

FOR I : =1 TO K DO 
BEGIN 

SUM: =0. O; 

NKK:=NIC 13; 

FOR J : =1 TO NKK DO 

SUM: =SUM-*-AC I , J 3 *EXP ( ( J-l ) *LN (TC IT! ) ) 
PROC II: =SUM*EXP ( — LC IltTCITI) ; 

END; 

FOR I : =1 TO K DO 
REL I : =REL I +PROC II; 

END; 

IF RELI >1.0 THEN 
RELC IT] : =1 . O 
ELSE 

REL C I T I : =REL I ; 

IF CHSIM=1 THEN 
RAND; 

END; 
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PROCEDURE YAXIS; 
BEGIN 



IF 


X 

II 

O 

N-' 


THEN 


PR I NTVAR: =LINE; 


IF 


( ( X=0) 


AND 


(Y=20>> AND (TCI 3=0) THEN PR I NTVAR 


IF 


( X=— 3) 


AND 


(Y MOD 


2=0) THEN PRINTVAR: =’ 0’ ; 


IF 


(X=-3) 


AND 


(Y=20> 


THEN PRINTVAR: =' 1 ’ ; 


IF 


( X=— 2) 


AND 


(Y MOD 


2=0) THEN PRINTVAR:=’ . ’ ; 


IF 


(X=-l) 


AND 


( (Y=20) OR (Y=0) > THEN PRINTVAR: =’ 


IF 


( x=-l ) 


AND 


(Y=1S) 


THEN PRINTVAR: =’ 9’ ; 


IF 


( X=— 1 ) 


AND 


(Y=16) 


THEN PRINTVAR: =’ 8’ ; 


IF 


<X=— 1 ) 


AND 


( Y= 1 4) 


THEN PRINTVAR: =’ 7’ ; 


IF 


( X=-l ) 


AND 


( Y= 12) 


THEN PRINTVAR: =’ 6’ ; 


IF 


( X=— 1 ) 


AND 


( Y= 1 0) 


THEN PR INTVAR: =’ 5’ ; 


IF 


(X=-l ) 


AND 


(Y=8) 


THEN PRINTVAR: =' 4’ ; 


IF 


( X=— 1 ) 


AND 


(Y=6) 


THEN PRINTVAR: =' 3' ; 


IF 


( X=— 1 ) 


AND 


(Y=4) 


THEN PRINTVAR: =’2' ; 


IF 


1 

II 

X 


AND 


( Y=2) 


THEN PRINTVAR: =’ I ' ; 



END; 

PROCEDURE GRAPHICS; 

BEGIN 

COUNT : =1 ; 

FOR Y: =20 DOWNTO 1 DO 
BEGIN 

FOR X : =-3 TO 76 DO 
BEGIN 

IF COUNT< =TMAX THEN 
BEGIN 

VAL:=R0UND(20*RELEC0UNT3) : 

IF (VAL=Y) AND (X=TCC0UNT3) THEN 
BEGIN 

PR I NT VAR : = ASTER I SK ; 

COUNT: =C0UNT+1 ; 

END 

ELSE 

PR I NT VAR : =BLANK ; 

END 

ELSE 

PRI NTVAR : =BLANK ; 

YAXIS; 

WRITE (PR I NTVAR) ; 

END; 

END; 

Y: =0; 

FOR X : =— 3 TO 76 DO 
BEGIN 

IF (X MOD 5 =0) THEN PR INTVAR: =L INE; 

IF COUNT < =TMAX THEN 
BEGIN 

VAL: =ROUND (20*REL [COUNT! > ; 

IF (VAL=Y> AND (X=TCC0UNT3) THEN 
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BEGIN 

PR I NT VAR : =ASTER I SK ; 

COUNT: =COUNT+l ; 

END 

ELSE 

IF (X MOD 5 <> 0) THEN 
PRINTVAR: =DASH; 

END 

ELSE 

IF (X MOD 500) THEN 
PR I NT VAR : =D ASH ; 

IF ( X=— 3) THEN PRINTVAR: = ’ O’ ; 

IF (X=-2) THEN PRINTVAR: =’ . ’ ; 

IF ( X = — 1 ) THEN PRINTVAR: =’ O’ ; 

WRITE (PRINTVAR) ; 

END; 

WRITE (’ O 5 10 15 20 25 30 35 40 45 

50 55 60’); 

WRITELN ( ’ TIME — >’ ) ; 

END; 



BEGIN 

INPUT; 

OUTPUT; 

FOR IT : =1 TO TMAX DO 
BEGIN 

INITIALIZE: 

REPEAT; 

ONEATATIME; 

UNTIL I E=K+ 1 ; 

SU92; 

IF CHS I M= 1 THEN 

WRITELN ( ’ RELIABILITY ( T >’ , TC ITU : 3, ’ ) = ’ , REL C I T D : 8 : & , 

’ (SIM = ’ , RELSIMC IT D : 5: 3, ’ ) ’ ) 

ELSE 

WRITELN ( * RELIABILITY ( T > ’ , T E I T 1 : 3, 7 ) = ' , RELE IT3 : S: 6) 

END; 

WRITELN; WRITELN; 

ATTENTION; 

WRITELN (’DO YOU ALSO WANT A GRAPHICAL PRESENTATION OF 
THE DATA?’ ) ; 

WRITE (’HIT "l" FOR YES OR "0" FOR NO : ’); 

READLN ( CHGR ) ; 

WRITELN; 

WRITELN; 

IF CHGR= 1 THEN 
GRAPHICS; 

DECIDE2; 

END. 
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APPENDIX B 



PROGRAM MI XED_CONVOLUT I ON; 

USES TRANSCEND, APPLESTUFF; 

CONST 

ASTERISK^ *’ ; 

BLANK = 7 7 ; 

DASH =' — 7 ; 

LINE = 7 ! 7 ; 

MAX =10; 

TYPE 

ARR = PACKED ARRAY I 1 . . MAX I OF REAL; 

ARI = PACKED ARRAY! 1. .MAXI OF INTEGER; 

DARR = PACKED ARRAY C 1 . . MAX , 1 . . MAX I OF REAL; 

VAR 

AR; PACKED ARRAY: C. . MAX, O. . MAX , O. . MAXI OF INTEGER; 

A; PACKED ARRAYC 1 . . MAX , 1 . . MAX 3 OF REAL; 

SUMREL: PACKED ARRAY CO.. MAX I OF REAL; 

N; PACKED ARRAYCO. . MAXI OF INTEGER; 

NINIT,NI,T, 2A, INP: ARI ; 

PROB, RELTOT, LAMB; DARR; 

L, PRC, REL, P; ARR; 

I, J,K,BR, BRANCHES, COUNT, KK, IL, KKK, J J , IE; INTEGER; 

I T . 1 1 , FA , NN , NKK , NN I , CHGR , TMAX , Y, X . VAL ; I NTEGER ; 

SUM , REL I , LAM , DEN ; REAL ; 

PRINTVAR, CH: CHAR; 

(* L = ARRAY OF ALL LAMBDA *> 

<* PRO: 1 1 = PROBABILITY OF EACH BRANCH *) 

(» T = ARRAY OF ALL TIMES AT WHICH THE RELIABILITY WILL 
BE COMPUTED *) 

<* BRANCHES = NUMBER OF BRANCHES IN SYSTEM *> 

PROCEDURE INPUT; 

BEGIN 

WR I TELN ( 7 HOW MANY DISSIMILAR LAMBDAS ARE IN THE SYSTEM? 7 > 
READLN (K) ; 

WR I TELN ( 7 WRITE DOWN ALL THE LAMBDAS WITH AT LEAST ONE BLANK 
IN BETWEEN: 7 ) ; 

WRITELN( 7 (NUMBERS < 1 NEED A LEADING ZERO!) 7 ); 

FOR I : = 1 TO K DO 
READ (LI 1 1 ) ; 

WR I TELN (’HOW MANY TIME POINTS DO YOU WANT TO BE COMPUTED? 

( <= 20 ) : 7 ) ; 

READLN (TMAX) ; 

WR I TELN ( 7 AT WHAT TIME DO YOU WANT THE SURVIVAL FUNCTION TO 
BE EVALUATED? 7 ); 

WR I TELN ( 7 ( INTEGER TIMES BETWEEN 0-60 IN THE SAME UNITS AS 
THE LAMBDAS) 7 ); 

FOR I : =1 TO TMAX DO 
READ (TC II); 



42 



WRITELN; WRITELN; 

END; 

PROCEDURE ATTENTION; 

VAR 

DURATION, I : INTEGER; 

PITCH; ARRAY! 1. . 103 OF INTEGER; 

BEGIN 

PITCH! 13;= 15; 

PITCH! 23 : = 17; 

PITCH! 3 3 : = 19; 

PITCH!43:= 20; 

PITCH! 53 : = 22; 

DURATION; = 30; 

FOR I ; =1 TO 5 DO 

NOTE (PITCH! I 3, DURATION) ; 

END; 

PROCEDURE OUTPUT 1; 

BEGIN 

WRITE (’FOLLOWING ’,K:3,’ LAMBDAS WERE ENTERED 
FOR I : =1 TO K DO 
WRITE (L!I3:5:2) ; 

WRITELN; 

END; 

PROCEDURE 0UTPUT2 ; 

BEGIN 

FOR I : = 1 TO BRANCHES DO 
BEGIN 

FOR J ; = 1 TO MAX DO 
BEGIN 

FOR K: =1 TO MAX DO 
BEGIN 

IF (ARII, J,K3=0> AND (KOI) THEN 
ARI I, J,K3:=-1; 

END; 

END; 

END; 

END; 

PROCEDURE GETDATA; 

BEGIN 

FOR I : =0 TO MAX DO 
FOR J : =0 TO MAX DO 
FOR K:=0 TO MAX DO 
ARII, J,K3:=-1; 

NI03 : = 1 ; 

WRITELN; 

WRITELN (’HOW MANY BRANCHES?’); 

READLN (BRANCHES) ; 

I ; =0; 
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REPEAT 

Is-I+1; 

WRITELN (1:3,’. BRANCH: ’>; 

J : =0; 

REPEAT 

J:=J+1; 

WRITELN (J: 3, ’ . ELEMENT: ’ >; 

K: =0; 

REPEAT 
K: =K+ 1 ; 

READ ( ZACK3 ) ; 

ARCI, J,K3:=ZACK3; 

NCJ3:=K; 

UNTIL ARC I , J , K3=0; 

until < (arci,j-i,ncj-i33=o) and <arci, j,ncj33=o> > 
AND (ZAC 13=0); 

UNTIL I=BRANCHES; 

0UTPUT2; 

END; 

PROCEDURE COMPPROB; 

BEGIN 

FOR I : =1 TO MAX DO 
FOR J : =1 TO MAX DO 
PROBCI, J3:=-l; 

I : = 1 : 

J:=2; 

WHILE I < =BRANCHES DO 
BEGIN 
DEN: =0; 

K: =2: 

COUNT: =1; 

WHILE ARC I , J , K3 >0 DO 
BEGIN 

WHILE AR C I , J , K 3 >0 DO 
BEGIN 

DEN: =DEN+LC ARC I , J , K3 3 ; 

K:=K+1; 

END; 

PROBC I , COUNT 3 : =LC ARC I , J , 1 3 3 /DEN; 

COUNT : =COUNT+ 1 ; 

J : =J+2; 

K: =2; 

DEN: =0; 

END; 

I : =1+1 ; 

J : =2; 

END; 

FOR I : =1 TO BRANCHES DO 
BEGIN 
J : = 1 ; 

PC I 3 : = 1 ; 
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WHILE PROBE I, J 1 >0 DO 
BEGIN 

PC I 3 : =PC I 2 *PROBC I , J 3 ; 

J : =J + 1; 

END; 

END; 

WRITELN; 

END; 

PROCEDURE COMPLAMB; 

BEGIN 

FOR I : =1 TO MAX DO 
FOR J : =1 TO MAX DO 
LAMB C I , J 3 : =— 1 ; 

I : =1 ; 

J : = 1 ; 

WHILE I<=BRANCHES DO 
BEGIN 
LAM: =0; 

K:=l; 

COUNT : =1 ; 

WHILE ARC I , J , K3 >=0 DO 
BEGIN 

WHILE ARE I , J , K3 >0 DO 
BEGIN 

LAM: =LAM+LE ARC I , J , K3 1 
K : =K+ 1 ; 

END; 

LAMB C I , COUNT I : =LAM; 

COUNT :=C0UNT+1; 




LAM: =0; 

END; 

1 : = 1 + 1 ; 

J:=l; 

END; 

END; 

PROCEDURE SORT; 

VAR 

TEMP: REAL; 

SORTED: BOOLEAN; 

BEGIN 
REPEAT 
I : =1 ; 

SORTED: =TRUE; 

WHILE LE I + 1 D >0 DO 
BEGIN 

IF LCI1>LCI+13 THEN 
BEGIN 

TEMP: =LC 13; 
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LE ID: =L E I + 1 3 ; 

LE I +1 3 : =TEMP; 

SORTED: =FALSE; 

END; 

I:=I+1; 

END; 

UNTIL SORTED=TRUE ; 

LE I +1 D : =— 1 ; 

LE 1 + 23 : =— 1 ; 

END; 

PROCEDURE INITIALIZE; 

BEGIN 
KK : = 1 ; 

WHILE KK< =K DO 
BEGIN 

AEKK, 1 D: =1 . 0; 

NI EKKD : =1 ; 

J J : =1 ; 

WHILE JJ<=K DO 
BEGIN 

IF JJOKK THEN 
BEGIN 

AEKK. 1 3 : =AEKK, 13*LCJJ3/ (LE J J 3-LEKK3 ) ; 

END; 

JJ : =J J + l : 

END; 

KK:=KK+1 ; 

END; 

KK:=K; 

IE: =1 ; 

END; 

PROCEDURE SUMMATION; 

BEGIN 

SUM: =0. O; 

FOR I : =1 TO K DO 
BEGIN 

IF IOIE THEN 
BEGIN 

SUM: =SUM+AE I , 13HLE I 3/ (LE I3-LE IE3) ; 

IF NI E I 3 >=2 THEN 
BEGIN 
FA: =1 ; 

NKK : =NI E I 3 ; 

FOR II: =2 TO NKK DO 
BEGIN 

FA: =FAt ( 1 1 — 1 > ; 

IF L E I 3 >L E I E 3 THEN 

SUM: =SUM+AE 1 , 1 1 3 *FA*LE IE3/EXP ( 1 1 *LN ( 
LEI3-LEIE3) > 

ELSE 
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BEGIN 

IF (II MOD 2>=0 THEN 

SUM: =SUM+AC I , I I3*FA*LC IE3/EXP ( 

1 1 *LN ( ABS (LC I 3— LC IE3 ) ) ) 

ELSE 

SUM: =SUM+AC I , I I 3*FA*LC IE3/ (-EXP ( 
II*LN(ABS(LCI3-LCIE3> ) ) ) 

END; 

END; 

END; 

END; 

END; 

AC IE, 1 3:=SUM+AC IE, 1 3; 

END; 

PROCEDURE SUB1 ; 

BEGIN 
I : = 1 ; 

WHILE I<=K DO 
BEGIN 

IF I < > I E THEN 
BEGIN 

AC I, NIC I 33: =ACI ,NI CI33SLCIE3/ (LCIE3-LCI3); 

IF NIC I 301 THEN 
BEGIN 
J : =1 ; 

NNI : =NI C I 3— J ; 

WHILE J< = (NI C I 3 — 1 ) DO 
BEGIN 

ACI,NNI3:=(LCIE3*ACI,NNI3-NNI*ACI,NNI+13)/ 
(LC IE3-LC 13); 

J: =J+1 ; 

NNI : =NI C I 3-J ; 

END; 

END; 

END: 

1 : = 1 + 1 ; 

END; 

END; 

PROCEDURE ONEATATIME; 

BEGIN 

IF NINITC IE3=NI C IE3 THEN 
IE: =IE+1 
ELSE 
BEGIN 

N I C I E 3 : =N ICIE3 + 1; 

J : =0; 

REPEAT 

NN: =NI C IE3-J ; 

AC IE, NN3 : =LC IE3<ACIE, NN— 1 3/ (NN-1) ; 

J : = J + 1 ; 
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UNTIL NN=2; 

SUMMATION; 

SUB1; 

END; 

END; 

PROCEDURE SUB2; 

BEGIN 

RELI : =0. O; 

IF TCIT3=0 THEN 
RELI : =1 . O 
ELSE 
BEGIN 

FOR I : =1 TO K DO 
BEGIN 

SUM: =0.0; 

NKK : =NI C I 3 ; 

FOR J : = 1 TO NKK DO 

SUM: =SUM-*-AC I , J3*EXP ( (J-l) *LN(TCIT3) ) 
PROCI3:=SUM*EXP(-LCI3*TCIT3) ; 

END; 

FOR I : =1 TO K DO 
RELI : =RELI +PROC ID; 

END; 

IF RELI >1.0 THEN 
RELCIT3:=1.0 
ELSE 

REL C I T 3 : =REL I ; 

END; 

PROCEDURE PREPFORLOOP; 

BEGIN 

FOR I : =1 TO MAX DO 
NINITC I 3 : =-l ; 

KKK : = 1 ; 

I : = 1 ; 

IF L C 2 3 < O THEN 
BEGIN 

NINITC 1 3:=1 ; 

K: = 1 ; 

END; 

WHILE LC 1 + 1 3 >0 DO 
BEGIN 

IF LCI3<LCI+13 THEN 
BEGIN 

NINITC KKK 3 : = 1 ; 

K : =KKK ; 

KKK: =KKK+1 ; 

I:=I+1; 

IF LC 1+1 3<0 THEN 
NINIT CKKK3 : = 1 ; 

K: =KKK ; 
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IF 


( X=— 1 ) 


AND 


CY=8) 


THEN 


PRINTVAR 


= ’ 4 ’ 




IF 


( X =— 1 ) 


AND 


( Y=6 ) 


THEN 


PRINTVAR 


= ’3’ 


? 


IF 


( X=— 1 ) 


AND 


( Y=4 ) 


THEN 


PRINTVAR 


= ’2’ 


5 


IF 


( X=— 1 ) 


AND 


( Y=2) 


THEN 


PRINTVAR 


= ’ 1’ 


5 



END; 

PROCEDURE GRAPHICS; 

BEGIN 

COUNT : = 1 ; 

FOR V : =20 DOWNTO 1 DO 
BEGIN 

FOR X : =-3 TO 76 DO 
BEGIN 

IF COUNT < =TMAX THEN 
BEGIN 

VAL:=R0UND(20*RELCC0UNTI) ; 

IF (VAL=Y) AND < X=TCC0UNT3 > THEN 
BEGIN 

PR I NTVAR : = ASTER I SK ; 

COUNT ;=C0UNT+1 ; 

END 

ELSE 

PR I NTVAR : =BL ANK ; 

END 

ELSE 

PR I NTVAR : =BL ANK ; 

YAXIS; 

WRITE (PR I NTVAR) ; 

END; 

END; 

Y: =0; 

FOR X : =-3 TO 76 DO 
BEGIN 

IF (X MOD 5 =0) THEN PR I NTVAR: =LINE; 

IF COUNT < =TMAX THEN 
BEGIN 

VAL: =ROUND ( 20*REL C COUNT! ) ; 

IF (VAL=Y) AND ( X=T:C0UNT3 ) THEN 
BEGIN 

PRINTVAR: =ASTERISK; 

COUNT: =C0UNT+1; 

END 

ELSE 

IF (X MOD 500) THEN 
PRINTVAR: =DASH; 

END 

ELSE 

IF (X MOD 500) THEN 
PRINTVAR: =DASH; 

IF (X=-3) THEN PRINTVAR: =’ O' ; 

IF ( X=-2) THEN PRINTVAR: ; 

IF (X=-l) THEN PRINTVAR: =’ O’ ; 



49 



END 
ELSE 
BEGIN 
1 1 : = 1 ; 

COUNT : = 1 ; 

WHILE LE 1 1 3=L C 1 1 + 1 3 DO 
BEGIN 

COUNT: =COUNT+l : 

N I N I T E KKK 3 : =COUNT ; 

FOR IL: =1 1 TO 9 DO 
L E I L 3 : =L E I L + 1 3 ; 

LEMAX3: =-l ; 

K:=KKK; 

END; 

KKK: =KKK+1 ; 

I : = 1 1 + 1 ; 

IF (( (COUNT >1) AND (LC 1+1 3<0) AND <LEI3>0))) THEN 
BEGIN 

NINITEKKK3: =1 ; 

K: =KKK; 

END; 

END; 

END; 

END; 

PROCEDURE LOOP; 

BEGIN 

FOR I T : =1 TO TMAX DO 
BEGIN 

INITIALIZE; 

REPEAT 

ONEATATIME; 

UNTIL IE=K+1; 

SUB2; 

RELC I T3 : =RELE IT3 *PE BR3 ; 

RELTOTEBR. IT 3 : =RELE I T3 ; 

END; 

END; 

PROCEDURE VAX IS; 

BEGIN 



IF 


(X=0> 


THEN 


PR I NT VAR: =LINE; 


IF 


( (X=0> 


AND 


(Y=20) ) 


AND 


(TE 1 3=0) THEN PRINTVAR: =ASTERISK 


IF 


( X=— 3) 


AND 


(Y MOD 


2=0) THEN PRINTVAR: =’ O’ ; 


IF 


( X=— 3) 


AND 


(Y=20> 


THEN 


PRINTVAR: =’ 1 ’ ; 


IF 


( X=-2) 


AND 


(Y MOD 


2=0) 


THEN PR I NT VAR 


IF 


( X=-l ) 


AND 


( <Y=20) 


OR 


(Y=0) ) THEN PRINTVAR: = ’ 0’ ; 


IF 


(X=-l) 


AND 


(Y=1S) 


THEN 


PRINTVAR: =’ 9’ ; 


IF 


( X =— 1 ) 


AND 


<Y=16) 


THEN 


PRINTVAR: =’ 8’ ; 


IF 


( X=- 1 ) 


AND 


<Y=14) 


THEN 


PRINTVAR: =’7’ ; 


IF 


( X=-l ) 


AND 


( Y= 1 2 ) 


THEN 


PRINTVAR: =’6’ ; 


IF 


( X=— 1 ) 


AND 


(Y=10) 


THEN 


PRINTVAR: =’ 5’ ; 
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WRITE (PRINTVAR) ; 

END; 

WRITE <’ 0 5 10 15 20 25 30 35 40 45 

50 55 60’ > ; 

WRITELN ( ’ TIME — >’ ) ; 

END; 

PROCEDURE DECIDE; 

BEGIN 

WRITELN; WRITELN; 

WRITELN (’DO YOU WANT TO CONTINUE WITH MIXED CONVOLUTIONS? 
TYPE "U" : ’ ) ; 

WRITELN (’OR DO YOU WANT TO COMPUTE CONVOLUTIONS? 

TYPE "XP1" :’); 

END; 

BEGIN 

INPUT; 

OUTPUT 1 ; 

GETDATA; 

0UTPUT2; 

COMPPROB; 

COMPLAMB; 

FOR I : =1 TO MAX DO 
L 1 1 3 : =-l ; 

FOR BR: =1 TO BRANCHES DO 
BEGIN 
J : = 1 ; 

WHILE LAMB I BR , J 3 >=0 DO 
BEGIN 

IF LAMB! BR , 1 3 >0 THEN 
BEGIN 

LI J 3 : =LAMBCBR, J 3 ; 

J : =J+1 ; 

END 

ELSE 

BEGIN 

L I J 3 : =LAMB I BR , J + 1 3 ; 

J : =J+1 ; 

END; 

END; 

SORT; 

PREPFORLOOP; 

LOOP; 

END; 

FOR I ; = 1 TO TMAX DO 
BEGIN 

SUMREL 1 1 3 : =0 ; 

FOR BR: =1 TO BRANCHES DO 

SUMREL 1 1 3 : =RELTOT I BR , I 3 +SUMREL 1 1 3 ; 

WRITELN( ’RELIABILITY ( T > ’,TCI3:3,’ ) = 

SUMREL! 13:8:6) ; 
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END; 

WRITELN; 

ATTENTION; 

WRITELN (’DO YOU ALSO WANT A GRAPHICAL PRESENTATION 
OF THE DATA?’); 

WRITE (’HIT "l" FOR YES OR "O” FOR NO : ’); 

READLN ( CHGR ) ; 

WRITELN; 

WRITELN; 

IF CHGR= 1 THEN 
GRAPHICS; 

DECIDE; 

END. 
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